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Abstract 

We discuss the event rate in DeepCore array due to neutrino flux produced by annihilations and 
" decays of galactic dark matter. This event rate is calculated with a 10 GeV threshold energy, which 

q; is smaller than the threshold energy taken in previous works. Taking into account the background 

event rate due to the atmospheric neutrino flux, we evaluate the sensitivity of DeepCore array for 
probing dark matter annihilation cross section and decay time. The sensitivity studies include the 
annihilation modes XX ~^ bb, t + t~ , and uv, and decay modes x ~ > bb, t + t~ , u + u~, and 

vv. We compare our results with corresponding constraints derived from observations of WMAP, 
ACT and Fermi-LAT. 

PACS numbers: 14.60.Pq, 14.60.St 
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I. INTRODUCTION 



Many astrophysical observations indicate the existence of dark matter (DM). A good 
example of such observations is the measurement of rotation curves for stars and gas in 
spiral galaxies. On the other hand, the nature of DM remains to be unveiled. In this regard, 
many DM candidates have been proposed with the proposal of weak interacting massive 
particles (WIMPs) jl[ ^ the most popular among all candidates. The detections of DM 
shall test the idea of WIMPs. The detections of DM can be categorized into direct and 
indirect approaches. The former approach proceeds by observing the nucleus recoil as DM 
interacts with the target nuclei in the detector. The latter approach relies on detecting 
final state particles resulting from DM annihilations or decays. In this article, we focus on 
the indirect detection of DM through observing neutrinos produced by DM annihilations or 
decays in the galactic halo. 

The search for neutrinos coming from DM annihilations in galactic halo has been per- 
formed by IceCube sjj. Data obtained from IceCube 22-string configuration set the 90% 
CL. upper limit on DM annihilation cross section (av) ~ 10~ 22 cm 3 s _1 for XX ~~ y v channel 
at m x = 1 TeV 4[ , while the preliminary result of IceCube 40-string galactic-center analysis 
improves the above limit to 10 _23 cm 3 s~ 1 [5]. Based upon IceCube 40-string data selection 
for searching diffuse flux of astrophysical muon neutrinos 6j, a comparable constraint on 
the cross section of XX ~^ y v annihilations, which could take place in the core of Earth, is 
derived j3j for TeV range dark-matter masses. Furthermore, upper limits on the annihilation 
cross sections of XX ~~ >* and XX ~~ ^ t + t~ channels are also obtained. It is interesting 

to compare these upper limits with the required annihilation cross sections for the same 



channels for explaining the PAMELA data on positron fraction excess [8| and Fermi-LAT 
e + + e~ fluxes measurement As shown in Ref. the IceCube upper limit on XX ~^ t~ + t~ 
annihilation cross section is comparable to the required XX ~~ t + t~ annihilation cross sec- 
tion for explaining PAMELA and Fermi-LAT data. On the other hand, the IceCube upper 
limit on XX ~~ >" annihilation cross section is still too high to test the idea of using this 

mode to account for PAMELA and Fermi-LAT data. 

The IceCube sensitivity on DM search is expected to improve with the data from all 86 
strings analyzed. The analysis of DeepCore array data will further enhance the sensitivity. 
The DeepCore array 



10 



12( is located in the deep center region of IceCube detector. This 
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array consists of 8 densely instrumented strings plus the nearest standard IceCube strings. 
The installation of DeepCore array significantly improves the rejection of downward going 
atmospheric muons in IceCube and lowers the threshold energy for detecting muon track 
or cascade events to about 5 GeV. This muon rejection is crucial for IceCube to observe 
DM induced neutrino signature from galactic halo. In fact, it has been pointed out that 
the parameter range for XX ~~ >* channel preferred by PAMELA and Fermi-LAT data 



could be stringently constrained 



13Ml5|(see also discussions in Ref. (l6j) by the data from 



IceCube detector augmented with DeepCore array. In this work, we shall however focus 
on the low-mass DM instead of DM with its mass in the range preferred by PAMELA and 
Fermi-LAT. 

We note that previous analyses on DeepCore sensitivity [jj], [n| have set the threshold 
energy at (40 — 50) GeV for both track and cascade events. However, to take the full 
advantage of DeepCore array, it is desirable to include neutrino events in the energy range 
10 GeV < E u < 50 GeV. We have initiated such a study for track events 18[. In this work, 
we generalize the previous study to cascade events in IceCube DeepCore detector. Since we 
are only interested in low-mass DM, we only consider channels XX ~~ >* bb, t + t~ , fi + fi~, and 
vv for annihilations and channels x ~ ► bb, t + t~ and and vv for decays. The neutrino 

fluxes generated through DM annihilations or decays into tt, W + W~ and ZZ final states 
are not included in this analysis. 

This paper is organized as follows. In Sec. II we describe neutrino fluxes from DM 
annihilation/decay at the galactic halo for different halo profiles. In Sec. Ill we briefly 
describe our results on the background atmospheric neutrino fluxes taking into account 
neutrino oscillations. In Sec. IV we present our results on the projected five year sensitivity 
of DeepCore array on cascade events induced by DM annihilations and decays. In addition, 
we shall also compare our results with the up-to-date indirect detection limits from Fermi- 
LAT gamma ray data and cosmic microwave background (CMB) observations. Specifically, 
we shall compare the DeepCore sensitivities on DM annihilation cross section (av) and 
DM decay time with corresponding constraints obtained from gamma ray observations 19, 



2o|(see also analysis in Ref. 2l|) and those obtained from CMB anisotropy 22|, |23j based 
on one or both of the recent WMAP 7-year [24| and ACT 2008 '% 
conclusions are given in Sec. V. 



data. Discussions and 
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II. NEUTRINO FLUX FROM ANNIHILATIONS AND DECAYS OF DARK MAT- 
TER IN THE GALACTIC HALO 



The differential neutrino ffux from the galactic dark matter halo for neutrino flavor i can 



be written as 



26) 



d$„. Att (av) 



for the case of annihilating DM, and 
d$„ Afl 1 




R Q p% x J 2 (Afi) (1) 



dE Vi 4.n m x T x 




R qPq x Jx(Afi) (2) 



for the case of decaying DM, where R Q = 8.5 kpc is the distance from the galactic center 
(GC) to the solar system, p = 0.3 GeV/cm 3 is the DM density in the solar neighborhood, 
m x is the DM mass, r x is the DM decay time and dN^JdE is the neutrino spectrum per 
annihilation or decay for a given annihilation or decay channel F with a corresponding 
branching fraction Bp. For XX ~~ >* y v channel, we assume that two neutrinos are produced 
per annihilated DM pair and all neutrino flavors are equally populated. Thus the neutrino 
spectrum per flavor is a monochromatic line with dN v /dE = \8(E — m x ). On the other 
hand, the neutrino spectrum per flavor for \ — > vv channel is dN u /dE = %5(E — ^). The 
neutrino spectra dN^/dE for other channels are summarized in Refs 13, |27[. The quantity 



(av) is the thermally averaged annihilation cross section, which can be written as 

(av) = B(av) , (3) 



with a boost factor B 28|]. We set (av)o = 3 x 10~ 26 cm 3 s _1 , which is the typical annihilation 
cross section for the present dark matter abundance under the standard thermal relic scenario 
We treat the boost factor B as a phenomenological parameter. The dimensionless 
quantity J n (AQ) is the DM distribution integrated over the line-of-sight (l.o.s) and averaged 
over a solid angle Afl = 2n(l — cos^ max ), i.e., 

j.(m = ±[ dn/ ^(™y, ( 4, 

Ail J A n Am R Q \ P& J 



where p is the DM density at a specific location described by the coordinate (I, ip) with I 
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TABLE I: Parameters of the density profiles for DM halo. 



DM halo model R s in kpc p s in GeV/cm 3 



NFW 20.0 0.260 

Einasto 21.5 0.0538 

Isothermal 3.50 2.07 



the distance from the Earth to DM and ip the direction of DM viewed from the Earth with 
ijj = corresponding to the direction of GC. The distance r = \J R%, + I 2 — 2R ( 7 ] lco svl> is the 



:SUSY 



30J, Einasto 



distance from GC to DM. We computed the values of J n (Af2) with Dar 
galactic DM distribution, we consider Navarro-Frenk- White (NFW) 
and Isothermal profiles 34j. The functional forms of these profiles are given by 

R s f r N 
PnfwO) = Ps-p- ( 1 + 



29). 



or 



the 



31 



(5) 



PEm{r) = p s X exp 



a 



a 



0.17, 



(6) 



PlsolJJ 



Ps 



(7) 



\ + {r/R s r 

with values of R s and p s given in Table I. In all cases we impose the normalization p & = 
0.3 GeV/cm 3 , which is at r = 8.5 kpc. The comparison of different DM halo models is 
shown in Fig. 1. 

Neutrinos are significantly mixed through oscillations when they travel a vast distance 
across the galaxy. The observed flavor ratio of DM induced neutrinos is related to the flavor 



ratio at the source through the probability matrix P a p = P{vp — > v a ) |35H37j . In particular, 
the exact form of P a p in terms of mixing angles 6*^ and CP phase 5 is given in Ref. j^ . 
The mixing angles 823 and 612 has been well measured while newest results for # 13 from 



accelerator 



38 



and reactor experiments 39M4l| are also available. We take sin 2 623 = 0.396, 



sin 2 9 12 = 0.307, sin 2 #13 = 0.0241 and 5 = I.O871", which are best fit values of neutrino mixing 
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parameters from a recent global fitting 42] in the case of normal mass hierarchy. Therefore, 



$> 
$„ 



0.55$" +0.24$" +0.21$ 



0.24$^ +0.40$" + 0.36$,, 



0.21$ 



0.36$ 



o 

I'r ' 




0.43$ 



o 



(8) 



where $°. and are neutrino fluxes at the source and on the earth respectively. While 
the best fit neutrino mixing parameters differ slightly in the case of inverted mass hierarchy, 
they do not produce noticeable change in the above relation. 



III. ATMOSPHERIC NEUTRINO FLUXES 



We follow the approaches in 18|, |43|, |44J to compute the intrinsic atmospheric neutrino 
background fluxes. The flux arising from ir decays reads 

d 2 Nl(E,C,X) 



dEdX 



[En Q(E n 
dE N / dE n - 

E J E 



1— 7tt 



d v E w (l - j n ) 
I -t — Pir{E 7T ,X,X')—F Nn (ET T ,E Nj 

A N 



E m 



x exp 



X' 
A N 



{E 



N \ J -'N J i 



(9) 



where E is the neutrino energy, X is the slant depth in units of g/cm 2 , £ is the zenith 
angle in the direction of incident cosmic-ray nucleons, r n = m^/m 2 , d n is the pion decay 
length in units of g/cm 2 , A ^ is the nucleon interaction length while A^v is the corresponding 
nucleon attenuation length, and <Pn{En) is the primary cosmic-ray spectrum. We have 
4>n(E n ) = J2 a A(/)a(E n ) with A the atomic number of each nucleus. The spectrum of each 



cosmic-ray component is parametrized by 



45 



^ A (E N ) =Kx(E N + bexpl-c^/E^}y 



(10) 



in units of m 2 s x sr x GeV 1 . The fitting parameters a, K, b, c depend on the type of nucleus. 



They are tabulated in Ref. 46j. The function P 7r (£' 7r , X, X') is the probability that a charged 
pion produced at the slant depth X' survives to the depth X (> X') [^J. F^yf-EL, En) is the 



normalized inclusive cross section for iV + air — > ^ + Y , and is given in Ref. 43]. The kaon 



contribution to the atmospheric flux has the same form as Eq. (Q with an inclusion of the 
branching ratio B(K — > fiu) = 0.635 and appropriate replacements in kinematic factors as 




Fig. 1: Comparison of the dark matter density distribution, p(r), as a function of distance from 
the Galactic Center as described by the NFW, Einasto and Isothermal halo models. 

well as in the normalized inclusive cross section. The three-body muon decay contribution 
to the atmospheric flux is also included. The details are discussed in Ref. 44|. After 
summing the two-body and three-body decay contributions, we obtain the total intrinsic 
atmospheric muon neutrino flux. From the Fig. 1 of Ref. 18|, we note that angle-averaged 



atmospheric muon neutrino flux obtained by our calculation and that obtained by Honda et 



al. 



48| both agree well with AMANDA-II results [49] . 



The intrinsic atmospheric v T flux due to D s decays can be obtained by solving cascade 
equations 44|, |50|. One obtains 

d 2 N Ur (E,X) Z NDs Z DsUt exp(-X/A N )4> N (E N ] 



(11) 



where Z 



NN 



dEdX 1 - Z NN A N 

1 — Aat/Aat, Z NDa and Z DaVr are the Z moments defined in our previous 



work 18|. Finally, the atmospheric flux taking into account the neutrino oscillation effect 
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is given by 



dE 



dX 



d 2 N, 



Vt_ p 
, ^ ■ -I 1/ 



dEdX 



dEdX 



(12) 



where P Ua ^ U/3 is the i/ a — > i/g oscillation probability. Sub-leading contributions to atmo- 
spheric i/ M flux arising from z/^ — >■ z/ e and v e — > oscillations are not included in the above 
equation. We can write down the atmospheric v T flux in the similar way. 

IV. RESULTS 



In IceCube DeepCore, the track event rate for contained muons is given by 



dE„ 



dE u N A p ice V t r x 



while the cascade event rate is given by 

n Ejji ax /* -E'max 



dE„ 



da^(E^E M 
dE„ 



\v ->■ v) 



;i3) 



dE. 



L't h 



d^A^PiceK; 



da uN (E u ,E sh ) 
dE v ' ' 



dE, 



(14) 



where p; ce = 0.9 gem 3 is the density of ice, Xa = 6.022 x 10 23 g 1 is Avogadro's num- 
ber, V tr ~ 0.04 km 3 is the effective volume of IceCube DeepCore array for muon track 



events 
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151 ] and V^ asc ~ 0.02 km 3 is that for cascade events 15|, |5jJ, d& v /dE u is the 



neutrino flux arrived at IceCube, which is the sum of DM induced flux and the background 
atmospheric neutrino flux, E max is taken as m x for annihilation and rn x /2 for decay, E^ and 
Ef^ are the threshold energies for track events and cascade events respectively, dcr„$ / dE^ is 
the differential cross section of neutrino- nucleon charged-current scattering, and da V N / dE s h 
is the differential cross section for showers produced by neutrino-nucleon charged-current 
and neutral-current scatterings. In this work, we use differential cross sections da^/dE^ 
and d<j V N I dE s h given by Ref. 52j with CTEQ6 parton distribution functions. The atmo- 



spheric part of d§ Ue /dE Ue is taken from Ref. 



. We also set E 



th 



10 GeV. It should be 



noted that the value for V tr is an average effective volume based on the energy dependent Vt T 



discussed in Ref. 
Core detector. The updatec 
events are available in Ref. 



10J, while the value for V^ asc is just the instrumental volume of the Deep- 
effective volumes of DeepCore detector for track and cascade 
53]. While we adopt constant effective volumes for evaluating 
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10 2 



Fig. 2: The dotted line, dot-dashed line, solid line, and dashed line are the expected sensitivities 
on DM annihilation cross section by DeepCore detector for XX ~~ ^ bb, XX ~^ t + t~ , XX ~ > M > 
and XX ~^ v ^ channels, respectively. 

DeepCore sensitivities, we shall also estimate how much the updated effective volumes affect 
our results. 

As mentioned earlier, we consider neutrino ffuxes generated through the annihilation 
channels XX ~~ bb, t + t~ , and i>V, and the decay channels x ~~ bb, r + r~, yU + /i~ and 

vv. By computing the cascade event rates, we present in Fig. 2 the required DM annihilation 
cross section as a function of m x such that the neutrino signature from DM annihilations 
can be detected at the 2a significance in five years. In this plot, we take NFW profile for 
DM density distribution in the galactic halo, the shower threshold energy = 10 GeV 
and the cone half- angle VVax = 50°. Non-detection of DM neutrino signature would then 
exclude the parameter region above the curve at the 2a level. We have presented results 
corresponding to different annihilation channels. It is seen that the required annihilation 
cross section for 2a detection significance is smallest for XX ~^ y v channel and largest for 
the channel XX ~^ bb. 
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At this juncture, it is desirable to estimate effects of the updated effective volume on our 
sensitivity calculations. The energy dependent effective volume V casc (E) as given in Ref. [53 1 
is roughly 3 times smaller than the value 0.02 km 3 adopted in our calculation for E u = 10 
GeV. On the other hand, V casc (E) increases monotonically with energy with V casc (E) greater 
than 0.02 km 3 for E u > 40 GeV. For the annihilation process XX ~ > one has E v = m x . 
Hence the DeepCore sensitivity to this process should be better than that presented in 
Fig. 2, which starts from m x = 50 GeV. For XX ~~ >* A^A 4 " mode, m x = 50 GeV corresponds 
to E v ~ 20 GeV. At this energy, the value V^ asc = 0.02 km 3 overestimates the effective 
volume by about a factor of 2, thus the DeepCore sensitivity to (<r(xX ~~ f- L+ f 1 ~) v ) should 
be corrected by a factor of \/2 « 1.4. The correction factor gradually reduces to 1 as m x 
approaches to about 120 GeV, which corresponds to E v = 40 GeV. For m x > 120 GeV, the 
DeepCore sensitivity calculated with V casc (E) is better than that presented in Fig. 2. The 
correction factor for XX ~^ t + t~ is similar to that for XX ~~ >" The hadronic mode 

XX ~~ >* bb requires different correction factor for the same m x . However we shall not address 
such a correction here since DeepCore detector is relatively insensitive to XX ~~ ^ bb. 

Next, we show how the DeepCore sensitivity on DM annihilation cross section varies with 
the chosen cone half-angle and threshold energy for NFW DM density profile. Here we take 
XX ~^ channel for illustration. Figure 3 shows the required DM annihilation cross 

section (cr(xx ~^ f I+ f 1 ~) v ) f° r a 2cr detection in five years for different cone half-angle ?/> max . 
It is seen that the sensitivity is improved as ?/> max increases from 1° to 2° while it turns 
weaker as ^ max increases further. In the latter case, the signal increases slower than the 
background does. In the same figure, we also show the result for a higher threshold energy 
E^ = 50 GeV with a cone half-angle ^ max = 10° for comparison. This result is taken from 



Ref. [17( where f/Vaix = 10° is identified as the most optimal cone half-angle for constraining 
DM annihilation cross section at that threshold energy. For large m x , lowering i?* h from 
50 GeV to 10 GeV does not affect much the signal rate while increases significantly the 
atmospheric background event rate. Hence, the sensitivity on DM annihilation cross section 
becomes worse by choosing E^ = 10 GeV. One expects the situation turns opposite for 
m x approaching to the threshold energy. In fact, for m x < 100 GeV, one can see that the 
sensitivity on DM annihilation cross section obtained with £"* h = 10 GeV is always better 



than that obtained with E£ = 50 GeV. We note that the DeepCore sensitivities on other 



annihilation channels have similar cone half-angle and threshold energy dependencies. 
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Fig. 3: The required DM annihilation cross section (xx ~ * as a function of m x such that the 

cascade events induced by neutrinos from DM annihilations can be detected at the 2a significance 
in five years. Results corresponding to different t/W x are presented. For comparison, we also show 



the result with E^} 1 



50 GeV and 



10 c 



17|. 



After discussing how the DeepCore sensitivities on DM annihilation cross section vary 
with the chosen cone half-angle and threshold energy for NFW profile, we study the variation 
of DeepCore sensitivities on XX ~>* with ^ max for different DM density profiles. First, 

we present in Fig. 4 our expected DeepCore sensitivities corresponding to different ^ max 
for Einasto DM density profile. Similar to NFW profile, Einasto DM density distribution 
also has a cusp in the central DM region. We refer to this class of DM density profiles as 
the cusped profile. Therefore the DeepCore sensitivity on (cr(xx — > fi + fi~)v) with Einasto 
profile becomes poorer as VWx increases. In addition, we can see from Fig. 1 that the DM 
density of Einasto profile is higher than that of NFW profile between 0.1 kpc and 5 kpc 
from the GC, and both density profiles are almost identical beyond this range of distances. 
Hence the DeepCore sensitivity on (<j(xX f i+ f I ~) v ) with Einasto profile is better than that 
with NFW profile for = 50°. Next, we also present in Fig. 4 the expected DeepCore 



11 



10 s 



— 10" 



E 

o 



x 

CO 



A 



10 1 



Einasto, 1 
' Einasto, 50° 



• Einasto, 90° 

Isothermal, 1° 

Isothermal, 50° 

Isothermal, 90° 

NFW, 50° 




10 4 



m z [GeV] 



10 J 



Fig. 4: The expected DeepCore sensitivities corresponding to different ^ max for Einasto (3ll. l32| 



, Isothermal 



3Jand NFW 



301 ] DM density profiles. 



sensitivities for different ip max with Isothermal DM density profile. Because there is a core in 
the central DM density distribution of Isothermal profile, we refer to this class of profile as the 
cored profile. We note that the sensitivity to DM annihilation cross section with Isothermal 
profile improves as ip max increases from 1° to 50° while the sensitivity with Einasto profile 
behaves oppositely. This is because that the DM density distribution of Isothermal profile 
maintains flat for a much longer distance from the GC as compared to the cusped profile. 
However, as VWx increases further, the DeepCore sensitivity to DM annihilation cross section 
becomes poorer even for Isothermal profile since the factor J 2 (Af2)Af2 is proportional to 
the square of DM density. Finally, DeepCore sensitivities to DM annihilation cross section 
with Isothermal profile are poorer than those with cusped profiles for the same ip max . This 
results from the fact that DM densities of cusped profiles around GC are much larger than 
that of Isothermal profile in the same region. 

In Fig. 5, we compare DeepCore sensitivities to DM annihilation cross section through 
measuring cascade events and track events. We note that is the dominant flavor of atmo- 
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Fig. 5: The required DM annihilation cross section (xx ~~ as a function of m x such that 
the neutrino signature from DM annihilations can be detected at the 2a significance in five years 
for track and cascade events. 

spheric neutrinos above few tens of GeV and the neutrino-nucleon cross section are almost 
the same for all flavors. Therefore, comparing with the result of track events, the signal 
to background ratio is enhanced in cascade events because only produces the cascade 
events through the neutral-current interaction, which is lower in cross section than that of 
charged-current interaction. Hence cascade events in general provide better sensitivities to 
DM annihilation cross section than those provided by track events. However, we can see 
in Fig. 5 that the sensitivity curve of track events for ip m3X = 2° is close to that of cas- 
cade events. In other words, for DM annihilations, the sensitivity provided by track events 
with sufficiently small cone half-angle can be comparable to the sensitivity provided cascade 
events. 

It is interesting and essential to compare our results with constraints obtained from 



gamma-ray astronomy and cosmology. Fermi large area telescope (Fermi-LAT) 54] is a 



pair-conversion telescope that explores the gamma-ray sky in the 20 MeV to 300 GeV range 
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Fig. 6: Cross section limits on DM annihilation channel \X ~~ M A* ■ The dot-dashed line is 
the lcr upper bound on annihilation cross section for low-mass WIMP [22|] obtained by WMAP7 



data. The dot-dot-dashed line is the CMB constraint obtained by using WMAP7 + ACT data at 



95% C.L. 



20J. The dotted line is the 



231 ] . The dashed line is the dSphs constraint at 95% C.L. 
constraint due to Fermi-LAT observations on the region \b\ > 10° plus a 20° x 20° square region 
centered at GC, assuming the NFW profile [2]. The solid line is our expected 2a sensitivity with 

ix — 50 . 



= 10 GeV and V: 



with unprecedented sensitivity. In a recent work, Fermi-LAT collaboration derive constraints 
on WIMP annihilation or decay into various final states which produce a continuous photon 
spectrum [if]]. These constraints are based upon the measured inclusive photon intensity 
spectrum from 4.8 GeV to 264 GeV obtained from two years of Fermi-LAT data over the 
region \b\ > 10° plus a 20° x 20° square region centered at GC with point sources removed. In 
Fig. 6, the dotted line is cross section upper limit on DM annihilation channel xx — y p + fi~ 
from the diffuse gamma-ray spectrum with NFW profile. It is taken from Ref. 19|] with a 
rescaling factor (4/3) 2 applied since we have adopted a local density of p & = 0.3 GeV/cm 3 
while Fermi-LAT analysis uses p & = 0.4 GeV/cm 3 . Our expected 2a sensitivity on XX ~^ 
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Fig. 7: Cross section limits on the DM annihilation channel XX ~^ bb. The dashed line is the dSphs 
constraint at 95% C.L. [20(. The dotted line is the constraint due to Fermi-LAT observations for 



the region \b\ > 10° plus a 20° x 20° square region centered at the GC, assuming the NFW profile 



19j. The solid line is our expected 2a sensitivity with Ef^ = 10 GeV and ip 



50°. 



fi + fi~ annihilation cross section by the DeepCore detector with ip m3X = 50° is plotted for 
comparison. We can see that our expected 2a sensitivity with T/> max = 50° is slightly stronger 
than this Fermi-LAT constraint, even if we calculate the sensitivity with V casc (E) given by 



Rcf. 



53]. 



Constraints on DM annihilation cross section were also obtained from cosmology and 
gamma-ray observations on dwarf spheroidal satellite galaxies (dSphs) of the Milky Way. En- 
ergy injection from DM annihilation at redshift 100 ^ z ^ 1000 affects the CMB anisotropy. 
This is because the injected energy can ionize the thermal gas and modify the standard 
recombination history of the universe. The updated CMB constraints on DM annihilation 



cross sections are derived in Refs. 



22] and 23| , where the former work is based on the recent 



WMAP 7-year data [2j while the latter one combines WMAP 7-year and ACT 2008 [25] 
data. We present the above two CMB constraints on XX ~* annihilation cross sec- 
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Fig. 8: The dotted line, the dot-dashed line, the solid line, and dashed line are the expected 
DeepCore sensitivities to DM decay time through cascade events from decay channels % — > bb, 
X — > t + t~ , x ~^ /^ + ^~) an d X ~^ vv, respectively. 

tion in Fig. 6. One can see that the solid line is higher than the dot-dot-dashed line (95% 
C.L.) by roughly a factor of 2. Thus our expected sensitivity is slightly weaker than the 
CMB constraint. However, it should be noted that the DM annihilation cross section could 
be velocity dependent. Hence, a model dependent extrapolation on DM annihilation cross 
section might be required to compare the constraint on (av) at redshift 100 ^ z < 1000 to 
that at the present day universe 55]. 

Dwarf spheroidal satellite galaxies of the Milky Way are DM-dominated systems which 
do not have active star formation or detected gas content 



56 



57] . Satellite galaxies are 



among the best targets to search for DM signals in gamma rays because of small background 
from astrophysical sources and a favorable signal to noise ratio. In Ref. {20)], Fermi-LAT 
collaboration derive upper limits on DM annihilation cross sections by applying a joint 
likelihood analysis to 24 months of data from 10 satellite galaxies with uncertainties on 
the dark matter distributions in the satellite galaxies taken into account. We present the 
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Fig. 9: The required DM decay time (x —> ^ + as a function of m x such that the neutrino 
signature from DM decays can be detected at the 2a significance in five years for cascade events. 
Results corresponding to different Vma? are presented. For comparison, we also show the result 



with Efi = 50 GeV and ip max = 50° 



121- 



dSphs constraint on XX ~~ >" ^ + ^~ annihilation cross section at 95% C.L. in Fig. 6. For m x > 
300 GeV, the expected DeepCore sensitivity curve is below that set by dSphs constraint. 
However, the dSphs constraint is stronger than the expected DeepCore sensitivity for m x < 
300 GeV, and even stronger than CMB constraint for m x < 100 GeV. 

We also present Fermi-LAT's dSphs constraint and galactic gamma ray constraint on 
XX ~~ bb annihilation cross section in Fig. 7. By comparing Figs. 6 and 7, we note that 
Fermi-LAT data gives more stringent constraint on XX ~^ bb mode than its constraint on 
XX ~^ mode. This is in contrast to the DeepCore case as one can see from Fig. 2. 

While the neutrino spectrum through XX is harder than that through XX ~^ bb, 

the gamma ray spectra through the above annihilations behave differently. The gamma ray 
spectrum from XX ~^ bb dominates over that from XX ~^ for most of the range of 

x = Ej/m x 58]. The expected 2a sensitivity on XX bb channel by the DeepCore detector 
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Fig. 10: The required DM decay time {x —> as a function of m x such that the neutrino 

signature from DM decays can be detected at the 2a significance in five years for track and cascade 
events. 

with ?/> max = 50° is also shown in Fig. 7. One can see that the expected DeepCore sensitivity 
on this channel is weaker than all existing constraints presented here. 

In addition to studying DeepCore sensitivities on DM annihilation channels, we also study 
sensitivities on DM decay time for \ ~~ >* bb, t + t~ , and v~P channels. Figure 8 shows 

the required DM decay time for reaching 2a detection significance in five years on cascade 
events from each channel with the threshold energy E^ = 10 GeV and a cone half-angle 
^max = 50°. Non-detection of such a signature would then exclude the parameter region 
below the curve at the 2a level. We can see that x — )• z/77 channel requires the lowest decay 
rate (longest decay time) to reach the 2a detection significance in five years of DeepCore 



data taking. We note that the energy-dependent effective volume V casc (E) [53[ downgrades 
the sensitivity to x ~ * y v decay time by roughly a factor of 1.5 at m x = 30 GeV, which 
corresponds to E v = 15 GeV. However, the sensitivity to \ ~ >" y u obtained with V casc (E) is 
better than that presented in Fig. 8 for m x > 80 GeV, which corresponds to E v > 40 GeV. 
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Fig. 11: The dot-dashed line and dashed line are the decay time constraints for x ~ > 
and x ~ * bb channels due to Fermi observations for the region |6| > 10° plus a 20° x 20° at GC, 
assuming the NFW profile Q. The solid line and dotted line are our expected 2a sensitivities 
for x — * an d X bb channels with = 10 GeV and ip ma , x = 50°. 

For x ~~ mode, the sensitivity at m x = 30 GeV is lower by a factor of \/3 ~ 1.7 by 
applying V C3SC (E), since such a m x corresponds to lowest possible neutrino energy E v = 10 
GeV. However, V casc (E) shall enhance the sensitivity for m x > 240 GeV, which corresponds 
to E v = 40 GeV, assuming E v ~ m x /6. The correction factor for x ~^ t + t~ is similar to 
that for x ~~ while the hadronic mode x ~ * bb requires different correction factor for 
the same m x . Once more, we do not address such a correction as the DeepCore detector is 
relatively insensitive to x ~ * bb- 

Next, we show how the DeepCore sensitivity on DM decay time varies with the chosen 
cone half-angle and threshold energy. We use x ~~ >" channel to illustrate these effects. In 

Fig. 9, we present the required DM decay time (% — > /i + yU~) as a function of m x for reaching 
2<j detection significance in five years for different cone half-angle ip max . The sensitivity curve 
rises as ip max increases from 1° to 50°. As ^ max increases in this cone half-angle range, the 
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DM event rate increases faster than that of atmospheric background. However, the curve 
falls with an increasing ip max for VWx > 50°, where the signal increases slower than the 
background does. We also show the required DM decay time for a 2a detection significance 
in five years with = 50 GeV and ?/Vax = 50° for comparison. It has been pointed 
out in Ref. 17j that ip max = 50° gives rise to the highest sensitivity on DM decay time for 



= 50 GeV. We note that the sensitivity on DM decay time is improved by lowering E^ 
from 50 GeV to 10 GeV for m x < 200 GeV. 

It is of interest to compare sensitivities on DM decay time given by cascade events and 
track events. In Fig. 10, one can see that cascade events provide better sensitivity on DM 
decay time than that given by track events for the same E^ 1 and ^ max . In the same figure 
we also show the sensitivity given by track events for ^ max = 90°. Such a ?/> max renders the 
best sensitivity for track events. However this sensitivity is still poorer than those given by 
cascade events. 

In Fig. 11, we present decay time lower limits for \ — > /i + fi~ and \ — > bb channels obtained 
from the diffuse gamma-ray spectrum for the region \b\ > 10° plus a 20° x 20° square region 
centered at GC, assuming the NFW profile. They are taken from Ref. [lj3], with a rescaling 
factor 3/4 applied since we have adopted a local density of p = 0.3 GeV/cm 3 while Fermi- 
LAT analysis uses p Q = 0.4 GeV/cm 3 . In the same figure, we also show the expected 2a 
sensitivities to x — > an d X ~ > bb decay time by the DeepCore detector with ^ max = 50° 

for comparisons. For the DeepCore detector, the sensitivity to the x ~^ decay time 

is better than that to the x ^ bb decay time, since the neutrino spectrum in the former 
channel is harder than the one in the latter channel. On the other hand, the Fermi-LAT 
data in general gives more stringent constraint on x ~ ► bb decay time than its constraint 
on the decay time of x ~ This is because that the gamma ray spectrum from the 
former channel dominates over the one from the latter channel for most of the range of 
x = E 1 /m x , as the DM annihilation case. If DM decays predominantly into one can 

see that the decay time sensitivity expected at DeepCore is comparable to the constraint 
given by Fermi-LAT data. This conclusion is not altered by adopting the energy-dependent 
effective volume V casc (E) 5jJ. On the other hand, if DM decays predominantly to bb, the 
expected DeepCore sensitivity is much poorer than the constraint from Fermi-LAT data. 
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V. DISCUSSIONS AND CONCLUSIONS 



In this paper, we have evaluated sensitivities of IceCube DeepCore detector to neutrino 
cascade events induced by DM annihilations and decays in galactic halo. We focus on the 
scenario of small DM mass and the threshold energy for the cascade events is taken to be 
10 GeV. The event rate of background atmospheric neutrinos is calculated with — > v T 
oscillations taken into account for neutrino energies less than 40 GeV. The signal event 
rate is calculated by taking NFW profile for DM density distribution in the galactic halo. 
Among all DM annihilation and decay channels, the annihilation mode XX ~~ >* y v an d the 
decay mode x — z/// provide the best search sensitivity, while the search sensitivity provided 
by the annihilation mode XX ~~ bb and the decay mode x ~^ bb is the poorest. 

It is important to compare the expected sensitivities of DeepCore detector to DM anni- 
hilation cross section and decay time with the existing constraints on the same quantities 
obtained by cosmology and Fermi-LAT gamma-ray observations. It is seen that the Fermi- 
LAT constraints on XX ~~ >" bb and x ~ * bb are much stronger than the expected DeepCore 
sensitivities to the same channels. Hence, if DM predominantly annihilates or decays into 
bb, the DeepCore detector is not expected to observe neutrino signature induced by DM 
in the galactic halo. For leptonic final states, one can see from Fig. 6 that the expected 
DeepCore sensitivity to XX ~~ >* annihilation cross section is slightly stronger than the 

Fermi-LAT constraint on the same channel based upon gamma-ray data from the galactic 
halo. However, the former is slightly weaker than both the dSphs constraint and the con- 
straint obtained from WMAP and ACT results on CMB anisotropy. In Fig. 11, we also 
see that the expected DeepCore sensitivity to x ~ >* decay time is comparable to the 

Fermi-LAT constraint based upon gamma-ray data from the galactic halo. From the above 
comparisons, there are still slight possibilities to observe DM induced neutrino signature 
from the galactic halo provided DM annihilates or decays predominantly into leptons. 

It should be noted that Fermi-LAT and CMB data does not directly set limits on XX ~~ y v 
or x ~ t y v modes with monochromatic neutrinos. There exist models (59), 60] in which 
XX ~~ y u an d X ~ y v are dominant annihilation and decay modes, respectively. In the 
annihilation case 59] for example, the DM candidate can be the lightest right-handed (RH) 



sneutrino in a U(1)b-l extension of the minimal supersymmetric standard model. RH 
sneutrinos annihilate into a pair of RH neutrinos. Each of these RH neutrinos then decays 
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into ordinary left-handed (LH) neutrino and a neutral Higgs boson while the decay of RH 
neutrino into charged final states l ± h T is typically forbidden in such a model. In the case 
that the mass difference between RH sneutrinos and RH neutrinos is small, RH neutrinos 
are produced non-relativistically by DM annihilations. Hence LH neutrinos produced by 
the decays of RH neutrinos are approximately monochromatic with an energy around half 
of the DM mass. Hence the sensitivity of IceCube DeepCore to this type of models can be 
read off from the XX ~~ >* y v curve in Fig. 2 with the shift ((av),m x ) — > {(cv), 2m x ). 

We like to stress that neutrino signatures induced by DM in the galactic halo can be 
more easily detected by neutrino telescope in northern hemisphere. The effective area for 
track events in KM3NeT as a function of neutrino energy has been estimated 61] . With this 
effective area, we also calculate the sensitivity of KM3NeT to DM annihilation cross section 
(<r(xX ->• in the galactic halo. For = 10° and m x = 200 GeV, the KM3NeT 

2cr sensitivity to (&(xx ~~ i n 5 years is comparable to the CMB constraint on this 
channel 23[. However, for m x = 1 TeV, CMB constraint gives (<j{xx ~~ >* f i+ f I ~) v ) no greater 
than 10~ 23 cm 3 s~ 1 whereas the KM3NeT 2a sensitivity on the same annihilation channel 
can reach to 1.5 x 10 _24 cm 3 s _1 in 5 years. 

In conclusion, we have made detailed comparisons between IceCube DeepCore sensitiv- 
ities and other existing constraints on various DM annihilation and decay channels. The 
prospect for DeepCore detector to observe neutrino signature induced by DM in the galactic 
halo has been discussed. We have also mentioned the possible improvement by KM3NeT on 
this observation. 
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